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Abstract: We present a rigorous proof of an ordering transition for a two-component two-dimen- 
sional antiferromagnet with nearest and next-nearest neighbor interactions. The low-temperature 
phase contains two states distinguished by local order among columns or, respectively, rows. Over- 
all, there is no magnetic order in accord with the classic Mermin- Wagner theorem. The method 
of proof employs a rigorous version of "order by disorder," whereby a high degeneracy among the 
ground states is lifted according to the differences in their associated spin-wave spectra. 



1. INTRODUCTION 

1.1 Background. 

For two-dimensional spin systems, the celebrated Mermin- Wagner theorem [32, 34] (and its ex- 
tensions [11,26]) precludes the possibility of the spontaneous breaking of a continuous internal 
symmetry. However, this result does not prevent such models from exhibiting phase transitions. 
For example, in the usual XY-model there is a low-temperature phase, known as the Kosterlitz- 
Thouless phase [28], characterized by power-law decay of correlations and, of course, vanish- 
ing spontaneous magnetization [22, 31]. The existence and properties of this phase have been 
of seminal importance for the understanding of various low-dimensional physical phenomena, 
e.g., 2D superconductivity and superfluidity, 2D Josephson arrays, 2D melting, etc. It it widely 
believed that no such phase exists for 0(n)-models with n > 3 although rigorous arguments for 
(or against) this conjecture are lacking. 

Of course, among such models there are other pathways to phase transitions aside from at- 
tempting to break the continuous symmetry. One idea is to inject additional discrete symmetries 
into the model and observe the breaking of these "small" symmetries regardless of the (global) 
status of the "big" one. As an example, at each r £ Z d (where d > 2) let us place a pair (a r , ir r ) 
of n-component unit-length spins whose interaction is described by the Hamiltonian 

H = -J\ ^2 (oy • a r i + 7IY • 7v) - J 2 ^(oy • 7r r ) 2 , (1.1) 

(r,r') r 
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where (r,r') denotes a pair of nearest neighbors on Z d and J 1; J 2 > 0. Obviously, this model 
has 0(n) symmetry (rotating all spins) as well as a discrete Z2 symmetry (relative reflection 
between the cr's and the 7r's). It is not hard to show that at low temperatures, regardless of the 
global status of the it's and 7r's, there is coexistence between a phase where the it's and 7r's 
are locally aligned with one another and one where they are locally antialigned. (Note that this 
is based purely on energy considerations — the said alignments are the only minimizers of the 
second term in the Hamiltonian.) A model similar to the one defined above was analyzed in [37] 
where the corresponding conclusions were indeed established. We remark that these results hold 
even if d = 2 (and even if n > 2). 

Another "circumvention" is based on the adaptation of the large-entropy methods to systems 
which happen to have continuous symmetry. These are distinguished from the more commonly 
studied systems in and of the fact that there is no apparent order parameter signaling the ex- 
istence of a low-temperature phase. The key idea dates back to [12, 29] where some general 
principles were spelt out that guarantee a point of phase coexistence. Let us consider an attractive 
system where there is an energetically favored alignment which confines the spin configurations 
to a small portion of the spin space. Suppose that there are many other less favored alignments 
with an approximately homogeneous energy. Under these conditions, a first-order transition at 
some (intermediate) value of temperature is anticipated. This kind of transition was established 
for specific systems (including the g-state Potts model) in [12,29], see also [38]. The general 
philosophy can easily be adapted to spin systems with a continuous symmetry, e.g., as in [2,7, 8] 
where some related problems were discussed. 

To illustrate these matters let us consider an example from [2]. Here we have a two-component 
spin of length one at each site of 1? which we parametrize by an angular variable 9 r E (— tt, it]. 
Let V(x) denote the function which equals negative one if \x\ < e and zero otherwise, and let 



where, of course, the arguments of V are interpreted modulo 2ir. Then, at some parameter 
value J = J t obeying e Jl « y/e, coexistence occurs between a phase where nearly all neigh- 
boring spins are closely aligned and one where, locally, spins exhibit hardly any correlation. We 
reiterate that the use of n = 2 and d = 2 is not of crucial importance for proofs of statements 
along these lines. Indeed, in [15, 16], similar results have been established in much generality. 

In all of the above examples a moment's thought reveals that no violation of the Mermin- 
Wagner theorem occurs. Indeed, this theorem does not preclude a phase transition, it only pre- 
cludes a phase transition which is characterized by breaking of a (compact) continuous internal 
symmetry. 

1.2 Foreground. 

The purpose of this note is to underscore another route "around" the Mermin- Wagner theorem. 
The distinction here, compared to all of the abovementioned, is that it may take the reader two 
moments to realize that our results are also in accord with the Mermin-Wagner theorem. Not 
unrelated is the fact that in our example the mechanism for ordering is relatively intricate. Let us 
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FIGURE 1. An example of the ground state of the Hamiltonian (1.3) on a finite grid. Here both 
sublattices exhibit Neel state with spins alternating between 30° and 210° on one sublattice and 
between 110° and 290° on the other. Any other ground state can be obtained by an independent 
rotation of all spins in each sublattice. 

go right to the (formal) Hamiltonian which reads 

r r 

Here r denotes a site in I? and the S r are unit-length two-component spins, i.e., S r G R 2 , 
with \S r \ = 1, for each r G I?. The vectors e x and e y are unit vectors in the x and y lattice 
directions while J (the overall interaction strength) and 7 (the relative strength of nearest neighbor 
couplings) are positive numbers. Notice the sign of the coupling — there is antiferromagnetism 
all around. 

In order to analyze the ground states, let us focus on the cases 7C 1, (Later we will only 
require 7 < 2.) Notice, especially in this limit, that the interaction splits the lattice into even and 
odd sublattices. For the ground-state problem, say in an even-sided finite volume with periodic 
boundary conditions, it is clear that both of the sublattices will be Neel (i.e., antiferromagneti- 
cally) ordered. However, once this Neel order is in place, it is clear that the energetics are insensi- 
tive to the relative orientation of the spins on the two sublattices. Specifically, the spin at any site r 
couples antiferromagnetically to the sum of S r+ g. x , S r+ £ y , S r _g. x and S r _£ y which, in any Neel 
state, is exactly zero. Thus we conclude that the set of ground states, i.e., the "order-parameter 
space," cf [33], of this model exhibits an 0(2) (g) 0(2) symmetry. 
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For convenience we will regard the first factor of O ( 2 ) <g> O ( 2 ) as acting on all spins and the sec- 
ond as acting on the relative orientations of (the spins on) the two sublattices. The upshot of this 
work (precise theorems will be stated in Section 2.1) is that, at small but positive temperatures, 
the order parameter space is reduced to Z 2 . Although the first 0(2) is restored as required by the 
Mermin-Wagner theorem, the remaining Z 2 is a remnant of the second 0(2). Consequently, at 
low temperatures, there are two Gibbs states: one where there is near alignment between nearest- 
neighbor spins in every lattice column and the other featuring a similar alignment in every lattice 
row. So the continuous 0(2) <g) 0(2) symmetry is evidently broken; we have Gibbs state in which 
all that acts is the single 0(2) factor. And all of this in two dimensions! 

Having arranged for the requisite two moments via procrastination, we will now reveal why 
this does not violate the Mermin-Wagner theorem. The answer is that the enhanced 0(2) (g> 
0(2) symmetry was never a symmetry of the Hamiltonian — this is both the hypothesis and the 
driving force of the derivations of the Mermin-Wagner theorem. Indeed, the large symmetry 
was only a symmetry of the ground state space and as such there is no a priori reason to expect 
its persistence at finite temperatures. So everything is all right. To further confuse matters, let 
us remark that although the "Z 2 remnant" — the one that does get broken — was not an internal 
symmetry of the Hamiltonian, it is, somehow, more organic than the 0(2) group that contained 
it. This particular Z 2 may be interpreted as the natural enactor of one of the lattice symmetries 
(here a 90° -rotation) which are typically associated with antiferromagnets. 

The last observation is supported by the fact that there is an order parameter associated with 
the above phase transition. Indeed, consider the object 

n r = (S r+6x -S r+6y )-S r (1.4) 

whose expectation is zero at sufficiently high temperatures and non-zero (in appropriate states) 
at low temperatures. (In another context, this sort of symmetry breaking has been referred to as 
Ising nematic ordering [1,27].) To summarize (in case all of this has been confusing), here we 
have a true long-range order but we avoid conflict with the Mermin-Wagner theorem because the 
0(2) (g) 0(2)-symmetry was never a true symmetry of the model. 

1.3 Order by disorder. 

In accordance with the title, the mechanism behind this ordering is called "order by disorder" (or, 
in the older vernacular, "ordering due to disorder"). This concept is, as of late, extremely preva- 
lent in the physics literature; most of the recent work concerns quantum large-,? systems where 
finite S plays the role of thermal fluctuations, but the origin of this technique can be traced to the 
study of classical systems, see [39,40] and [24]. In particular, in the latter reference, it is exactly 
the present model that was studied and this has since been referred to as the canonical model 
of order by disorder. The key words are "spin waves" and "stabilization by finite temperature 
excitations," neither of which should be unfamiliar to the mathematical physicist but which, until 
recently, have not been exploited in tandem. 

Let us proceed with the key ideas; we will attend to the obligatory citations later. For ease of 
exposition, let us imagine that somehow even at finite temperatures the two sublattices remain 
locked in their Neel states. Thus there is an angle, 0*, which measures the relative orientation 
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of the states on the two sublattices. Next we perform a spin-wave calculation to account for the 
thermal perturbations about the ground state with fixed cf>*. Although said instructions may have 
profound implications in other contexts, for present purposes this simply means "pitch out all 
interactions beyond quadratic order and perform the resulting Gaussian integral." The upshot of 
such a calculation is a quantity, the spin-wave free energy, which should then be minimized as a 
function of (f>*. As we will see this minimum occurs exactly when the states are either horizontally 
or vertically aligned, i.e., (f)* = 0° or <p* = 180°. 

The reader may question the moral grounds for the working assumption of finite temperature 
Neel order which is the apparent basis of the spin-wave calculation. Of course, the cheap way 
out — the final arbitrator — is the fact that herein is a rigorous proof. However, the spin-wave con- 
clusions are not so difficult to understand. Foremost, we reemphasize that the outcome is decided 
purely on the basis of free energetics. A cursory examination of the calculational mechanics 
then reveals that in fact only two ingredients are really needed. The first is that Neel order is 
present locally — which is certainly true at very low temperatures. The second boils down to the 
statement that the thermodynamic properties in these sorts of magnets are unaffected — to first 
approximation — if the system is restricted to configurations that have magnetic order. In particu- 
lar, the long wave-length excitations which are ultimately responsible for the break-up of ordering 
in two dimensions contribute insignificantly to the free energy. 

Now let us discuss the historical perspective of the present paper. The first phase in under- 
standing this sort of problems is coming to terms with the degeneracy of the ground-state space. 
When these situations arise, there is a selection at finite temperature according to the ability that 
each state has to harbor excitations. The simplest cases, namely a finite number of ground states 
and a small effective activity (e.g., a large "mass") for the excitations have been understood by 
physicists for a long time and are now the subject of essentially complete mathematical theo- 
rems [35,36,41]. Many interesting situations with infinitely many ground states were introduced 
in late 1970s and early 1980s, see, e.g., [17,40]. Here intricate and/or mysterious calculations 
are invoked to resolve the degeneracies — often resulting in phantasmagorical phase diagrams, see 
e.g. [18] — but the upshot in these situations is pretty much the same. In particular, with excruci- 
ating effort, some cases can now be proclaimed as theorems [6, 10]. However, the cornerstone of 
any systematic analysis (either mathematical or physical) is the existence of a substantial gap in 
the energy spectrum separating those excitations which resolve the ground-state degeneracy from 
the excitations that are readily available to all ground states. 

The degenerate ground-state problems look very different for the classical 0(n)-spin models. 
Indeed, the continuous nature of the spins in combinations with their internal degrees of freedom 
almost inevitably lead to a gapless excitation spectrum. Although this sounds a lot harder, the 
necessary computations turn out to be far more palatable. To our knowledge, the first such exam- 
ple, studied in [39], was a frustrated FCC antiferromagnet. The system is quite similar to the one 
discussed here but with the ordering caused, mostly, by quantum effects. In [24], studying exactly 
the model in (1.3), it was demonstrated that these techniques also apply to classical systems. In 
the present work we will transform these classical finite-temperature derivations into a mathemat- 
ical theorem. The proofs are quite tractable; all that is really required are some error estimates for 
the Gaussian approximations and a straightforward contour argument. To ease our way through 
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the latter we will employ the method of chessboard estimates. In some concurrent work [3,4], 
a similar analysis is used to resolve some controversies concerning models of transition-metal 
oxides. However, in these "TMO-problems," the ground-state spaces have additional intricacies 
so the beauty and simplicity of the method is obscured. 

To make our historical perspective complete let us also relate to the existing mathematical work 
on systems with continuous spins. A general approach to continuous spins with degeneracies has 
been developed in [13,42]. Here the method of resolution appears to be not terribly dissimilar to 
ours; e.g., there are quadratic approximations, Gaussian integrals, error estimates, etc. However, 
only a finite number of ground states are considered and we suspect that a detailed look at the 
"curvature conditions" will reveal that again there is a substantial mass gap in the excitation 
spectrum. Finally, from an earlier era, there are the methods based on infrared bounds [14, 19-21]. 
However, the reflection symmetries required to get these arguments started do not seem to hold in 
the system defined by (1.3). And even if they did, due to the infrared divergence, this would only 
provide misleading evidence — a la Mermin- Wagner — that the model under consideration has no 
phase transition. 



2. Main results 

2.1 Phase coexistence. 

To state our results on phase coexistence in the model under consideration, we will first recall the 
concept of infinite-volume Gibbs measures. We begin with finite-volume counterparts thereof, 
also known as Gibbs specifications. Let S = (S\,Sa<>) be a spin configuration where S\ 
and S ac denote the corresponding restrictions to A and A c , respectively. Let Ha{Sa, Sac) be 

the restriction of (1.3) to pairs of sites at least one of which is in A. Then we let [if 1 ^ be the 
measure on configurations in A defined by 

/4 SAc) (dSA) = Za(Sac) ^A(d5 A ). (2.1) 

Here S7a denotes the product Lebesgue measure on the unit circle, one for each r£ A. Following 
the "DLR-philosophy," see [23], the infinite- volume Gibbs measures are those measures on full 
configurations on 1? whose conditional probability in a finite volume A given the configuration 
in the complement is exactly the object in (2.1). 

In accord with the standard terminology, see [23], we will say that there is a phase coexistence 
for parameters J, 7 and (3 if there exists more than one infinite-volume Gibbs measure for the 
interaction (1.3) and inverse temperature [5. To adhere with mathematical-physics notation, we 
will refer to the Gibbs measures as Gibbs states and we will denote the expectations with respect 
to such states by symbol (— )p. 

Now we are in a position to state the main result of this paper. 

Theorem 2.1 Consider the model as defined above with fixed J G (0, 00) and 7 G (0, 2). Then 
there exists a (3o £ (0, 00) and a function (3 1— ► e(/3) satisfying e((3) — > as (3 — ► 00 such that 
the following holds: For each (3 > (3q there exist two distinct Gibbs states (— )^ and (— )^ 
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such that 




(2.2) 



(2.3) 



whenever r,r' G 71? are such thatr' = r + e Q . 

Let us informally describe the previous result. First, on both even and odd sublattice of Z 2 we 
have a (local) antiferromagnetic order. The distinction between the two states is that in (— 
the nearest-neighbor spins on I? are aligned in the x direction and antialigned in the y direction, 
while in (— )^ the two alignment directions are interchanged. In particular, it is clear that the 

order parameter rt r , defined in (1.4), has positive expectation in the z-state (—}^ and negative 

expectation in the y-state (—}^\ Since, as mentioned previously, Gibbsian uniqueness guarantees 
that (n r )j3 = at sufficiently high temperatures, we have a bone fide phase transition of the 
"usual" type. 

Despite the existence of multiple low-temperature Gibbs states, we emphasize that no claim 
has been made about the actual direction that the spins will be aligned to. On the contrary, we 
have the following easy corollary of the aforementioned Mermin- Wagner theorem: 

Theorem 2.2 Consider the model as defined above with J, 7 € M. fixed and let (— )p be any 
infinite-volume Gibbs state at inverse temperature (3. Then (— )p is invariant under the simulta- 
neous rotation of all spins and, in particular, {S r )p = Ofor allr G I?. 

The authors do not see any significant obstruction of Theorem 2.1 (appropriately modified) in 
the cases n > 2 and d > 2. For the case under consideration, namely, n = 2 and d = 2, it may 
be presumed that there is a slow decay of correlations at sufficiently low temperatures. Here it is 
conceivable that, with great effort, this could be proved on the basis of technology that is currently 
available [9,22, 30]. The anticipation is that for d > 3 and n > 2 there are actual sublattice Neel 
states while for d = 2 and n > 2 the decorrelations should be exponential. However, we do not 
expect to see a proof of any statement along these lines in the near future. 

2.2 Outline of the proof. 

We proceed by an informal outline of the proof of our main result (Theorem 2.1). The argument 
hinges on the following three observations: 

(1) Suppose A is a number that satisfies 



Then the (angular) difference of any typical pair of 7i<ytf-nearest neighbor spins will not devi- 
ate by more than A from the energetically optimal configuration. 
(2) In situations when (1) applies and under the additional assumption that A also satisfies 



f3JA 2 > 1. 



(2.4) 



/3JA 3 < 1, 



(2.5) 
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then all important contributions to the free energy of the system will come from a quadratic — 
or spin-wave — approximation to the Hamiltonian. 
(3) Finally, if F(rf>*) denotes the spin-wave free energy above the ground state where one sub- 
lattice is rotated by angle <p* relative to the other (see Fig. 1), then F ((/>*) is minimized only 

at <f>* = 0° or <j)* = 180°. 

(The mathematical statements corresponding to (1-3) above are formulated as Theorems 3.1 
and 3.2 in Section 3.1.) We observe that the necessary A as stipulated by (2.4-2.5) defines a 
running scale — not too big and not too small — which obviously tends to zero as j3 —> oo. 

Here is how these observations will be combined together to establish long-range order: We 
partition the lattice in blocks of side B. On the basis of (1) above, every block will with high 
probability exhibit a near ground-state configuration, which by (2-3) will have the sublattices 
either nearly aligned or nearly antialigned. Then we need to show that each of the two possibilities 
are stable throughout the entire system. For that we will resort to a standard Peierls' argument. 
Here the crucial observation (see Lemma 4.6) is that two "good" blocks with different type of 
alignment between sublattices are necessarily separated by a "surface" of "bad" blocks — that is 
those which either contain energetically charged pair of nearest-neighbor spins or whose spin- 
wave free energy exceeds the absolute minimum by a positive amount. 

Appealing to chessboard estimates (see Section 4.1), the probability of a particular "surface" 
can be factorized — as a bound — into the product over the constituting blocks. It turns out that the 
energetically frustrated "bad" blocks are suppressed once 

/3JA 2 >log£, (2.6) 

while the entropically frustrated blocks are suppressed once the excess spin-wave free energy 
times B 2 is sufficiently large. Under the conditions (2.5-2.6) and B » 1 the entropy of the 
above "surfaces" can be controlled. The desired phase coexistence then follows by standard 
arguments. 

A couple of remarks are in order: Due to the perfect scaling properties of Gaussian distribu- 
tions the suppression extracted from the spin-wave calculation is independent of (5 — the desired 
decay is achieved solely by choosing B sufficiently large. Large [5 is needed only to suppress 
large deviations away from the "perfect" ground states. Notwithstanding, for (technical) ease of 
exposition we will have to make B increase slowly with f}\ see (4.19) for the precise relation 
of A, B and 0. 

The various steps of the proof are laid out in the following order: In Section 3 we carry out 
the harmonic approximation and provide the needed control of the spin-wave free energy. In 
Section 4 we invoke chessboard estimates and some straightforward bounds to control the contour 
expansion. The actual proof of Theorem 2.1 comes in Section 4.3. 

3. Spin-wave calculations 



As mentioned above, the underpinning of our proof of the main result is (the outcome of) a 
spin-wave free-energy calculation. This calculation involves simply working with the harmonic 
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approximation of the Hamiltonian (1.3) for deviations away from a fixed ground state. The cal- 
culation itself is straightforward although special attention must be paid to the "zero mode." For 
reasons that will become clear in Section 4 — and also to make discrete Fourier transform readily 
available — all of the derivations in this section will be carried out on the lattice torus of L x L- 
sites. Here, for technical convenience, we will restrict L to multiples of four so that we can assure 
an equal status of the two Neel states. 

3.1 Harmonic approximation. 

We will begin by an explicit definition of the torus Hamiltonian. Here and henceforth we will 
parametrize the spins by angular variables = (9 r ) which are related to the S/s by the usual ex- 
pression S r = (cos 9 r , sin 9 r ). (Of course, the 9/s are always to be interpreted only modulo 2tt.) 
Up to irrelevant constants, the corresponding torus Hamiltonian Hl can then be written as 

H L (0) = jJ2{2 + COS(0, - e r+6x+6u ) + CO S (9 r - r+ g x -e y )} 
reT L 

reT L 

The spin-wave calculations are only meaningful in the situations where each of the sublattices is 
more or less aligned with a particular Neel state. To describe the overall and relative orientation 
of the spins on the even and odd sublattices we will need two angles 9* and (f>*, respectively. 
Depending on the parities of the coordinates of r, we will write the 9 r for r = (x, y) in terms of 
the deviation variables $ r as follows: 

9*, x,y-even, 

0* + (f>*, x-odd, y-even, ^ 

9* + it, x, y-odd, 

9* + 4>* + 7r, x-even, y-odd. 

Obviously, only the relative angle (f>* will appear in physically relevant quantities; the overall 
orientation 9* simply factors out from all forthcoming expressions. 

The principal object of interest in this section is the finite- volume free energy, which will play 
an important role in the estimates of "entropically-disfavored" block events in Section 4. For 
reasons that will become clear later, we will define this quantity by the formula 

F LA m = ~Iog / e-^V(«) © LV2 II (3-3) 

Here d9 r is the Lebesgue measure on unit circle and xl,a{0) = Xl,a(®'> ■> ^*) l& tne indicator 
that the deviation quantities #, defined from 6 as detailed in (3.2), satisfy | -d r \ < A for all r G T l. 
The factors of ^- have been added for later convenience. 

The goal of this section is to (approximately) evaluate the thermodynamic limit of the quan- 
tity -Fl, A (</>*) and characterize where it achieves its minima. As is standard in heuristic cal- 
culations of this sort, we will first replace the Hamiltonian (3.1) by its appropriate quadratic 
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approximation. We will express the resulting quantity directly in variables $ r : 

+ ^ 7 «*(**) £{(*," ^+ej 2 + (*r " ^ + e H ) 2 }- (3.4) 

ref L 

This approximation turns the integral in (3.3) into a Gaussian integral. As we will see later, here 
the indicator in (3.3) can be handled in terms of upper and lower bounds which allow "diago- 
nalization" of the covariance matrix by means of Fourier variables. The result, expressed in the 
limit L — > oo, is the following momentum integral: 

F W) = \I -^logZW), (3-5) 

where 

D k ((/)*) = |1 - e i(fcl+fc 2 )|2 + |j _ e i(fci-fc 2 )|2 +7COS (0*)(|! _ e ifcl|2 _ |J _ gifci|2^ (3 6) 

Here fci and /c 2 are the Cartesian components of vector k. The quantity F ((/>*) has the interpre- 
tation — justified via the preceding derivation — as the spin-wave free energy. As is checked by 
direct calculation, for 7 £ [0, 2) we have Dk((j)*) > almost surely with respect to the (normal- 
ized) Lebesgue measure on [—it, it] 2 . 

Having sketched the main strategy and defined the relevant quantities, we can now pass to the 
statements of (admittedly dry) mathematical theorems. First, we express the conditions under 
which the above approximate calculation can be performed: 

Theorem 3.1 Given e > and 7 6 [0, 2), there exists 8 = <5(e, 7) > such that if (5 J, A and 5 
satisfy the bounds 

PJA 3 < 5 and (3JA 2 > 1/5, (3.7) 

then 

limsup|F LiA (^) - F(<f>*)\ < e (3.8) 

L^oo 

holds for every (ft* £ (— 7r,7r]. 

The proof is postponed to Section 3.2. Having demonstrated the physical meaning of the 
function (ft* F(<fi*), we can now characterize its absolute minimizers: 

Theorem 3.2 For all 7 G (0, 2), the absolute minima of function (ft* ^ F ((/>*) occur (only) at 
the points (ft* = 0° and (ft* = 180°. 

Proof. The proof is an easy application of Jensen's inequality. Indeed, let a 6 [0, 1] be the number 
such that 2a — 1 = cos(^>*). Then we can write 

D k {<ft*) = aD k (0°) + (1 - a)£> t (180 o ). (3.9) 

Since -D&(0 o ) is not equal to D&(180°) almost surely with respect to dk (this is where we need 
that 7 > 0), the concavity of the logarithm and Jensen's inequality imply that F((f>*) > aF(0° ) + 
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(1 — a)F(180°) whenever 0,1. This shows that the only absolute minima that F can have 
are 0° and 180°. Now F is continuous (under the assumption that 7 < 2) and periodic, and so 
there exists at least one point in (— tt, tt] where it attains its absolute minimum. But F(0°) = 
F(180°) and so (f>* h-> F(4>*) is minimized by both </>* = 0° and ^ = 180°. □ 

3.2 Proof of Theorem 3.1. 

Throughout the proof we will fix J E (0, 00) and 7 E [0, 2) and suppress these from our notation 
whenever possible. Since everything is founded on harmonic approximation of the Hamiltonian, 
the starting point is some control of the error that this incurs: 

Lemma 3.3 There exists a constant c\ E (0, 00) such that the following holds: For any A E 
(0,oo), any 6>*,0* E (— 7r,7r] and any configuration 6 = (9 r ) of angle variables on Tt,, if the 
corresponding # = (i9 r ) satisfy \i} r \ < Afar all r E Tl, then 

\PH L (0) - < ci(1 + j)(3JA 3 L 2 . (3.10) 

Proof. We begin by noting that |# r | < A for allr E implies that |i? r — d r i \ < 2 A for all pairs 
of nearest and next-nearest neighbors r,r' E Tj,. This and the uniform bound 

Id 3 

cos(a + x) — (cos(a) — sin(a)x — \ cos(a)x 2 ) < , (3.11) 



show that, at the cost of an error as displayed in (3.10), we can replace all trigonometric factors 
in (3.1) by their second-order Taylor expansion in differences of "& r . Hence, we just need to show 
that these Taylor polynomials combine into the expression for 1l,4>*- 

It is easily checked that the zeroth order Taylor expansion in -& r exactly vanishes. This is a 
consequence of the fact that for $ = we are in a ground state where, as argued before, both 
sublattices can be independently rotated. This means we can suppose that 9* = <j>* = 0° in (3.1) 
at which point it is straightforward to verify that Hl{0) actually vanishes. Similarly easy it is to 
verify that the quadratic terms yield exactly the expression for Zl,<^*- It thus remains to prove that 
there are no linear terms in i? r 's. 

First we will note that all next-nearest neighbor terms in the Hamiltonian certainly have this 
property because there we have 9 r — 9 T +e; x ±t y ~ or 7r, at which points the derivative of the 
cosine vanishes. Hence we only need to focus on the nearest-neighbor part of the Hamiltonian — 
namely, the second sum in (3.1) — which we will temporarily denote by H™- Here we will simply 
calculate the derivative of with respect to i9 r : 

W HT{d) -=n= MOr+c* - Or) + Sin^+s, - Or) 



r 



- { sm(9 r - 9 r _ ix )+sm(9 r - 9 r _ iy ) } , (3.12) 

where the 9 r on the right-hand side should be set to the "ground-state" values. To make the 
discussion more explicit, suppose that r has both coordinates even. Then an inspection of (3.2) 
shows that the first sine is simply sm(<^>*) while the second sine evaluates to sin(0* + tt) = 
— sin(</>*). The net contribution of these two terms is thus zero. Similarly, the third and the fourth 
sine also cancel out. The other possibilities for r are handled analogously. □ 
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Using the harmonic approximation of the Hamiltonian, let us now consider the corresponding 
Gaussian equivalent of the integral in (3.3): 

reT L 

where d^ r is the Lebesgue measure on R and xl,a{'&) is the indicator that \& r \ < A for all 
r G Tl- Our next goal is to evaluate the effect of this indicator, which we will accomplish by 
proving an upper and lower bound on Ql,a(<^*)- We commence with the easier of the two, the 
upper bound: 

Lemma 3.4 For all (3 G (0, oo), all A G (0, oo) and all (p* G (— n, ir], 

limsup 1069 ^^^ -^). (3.14) 

Proof. The argument is relatively straightforward so we will be correspondingly brief. (A more 
verbose argument along these lines can be found in [3].) Pick a A > 0. We will invoke the 
exponential Chebyshev inequality in the form 

Xl,aW < e^ JXAL2 exp{-i/m £ K| 2 }. (3.15) 

r£T L 

Next we plug this bound into (3.13), diagonalize Zl,4>* by passing to the Fourier components 
= L^ 1 J2 r e7 L r &r& ir k and perform the Gaussian integrals with the result 

Here T£ = {27rL _1 (ni, 112) '■ ni = 1, 2, . . . , L} is the reciprocal lattice and ((f)*) is as defined 
in (3.6). The result now follows by taking logarithm, dividing by L? and invoking the limits L — > 
00 followed by A J, — with the last limit justified by the Monotone Convergence Theorem. □ 

The corresponding lower bound is then stated as follows: 

Lemma 3.5 For all ft G (0, 00), all A G (0, oo), all 0* G (— tt,tt] and all A > satisfying 
f3JA 2 X > 1, we have 



liminf 10 ^'^ 



>- W A) + log(l- * ), (3.17) 



l^oo L 2 ~ vr ' ; b V /3JA 2 A, 
w/iere F((f>*, A) is g/ve« fty j« integral as in (3.5) w/fft Dfr(<f)*) replaced by A + D^^cf)*). 

Proof. Again, we will be fairly succinct. Let A > 0. We begin by considering the Gaussian 
measure defined by 

r ^ m = 010b) -pf^w - e i"'i 2 }(f ) 12/2 n «* o-> 8 > 
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where Ql(,<P*,^) is an appropriate normalization constant. Let E\ denote the corresponding 
expectation. A simple bound shows that we have 

QlaW) > QU<t>*, A)E a (xa,l), (3.19) 

which reduces the desired estimates to two items: a calculation of the integral Ql(<I>*, A) and a 
lower bound on Ea(xa,l)- 

The first problem on the list is dispensed with similarly as in the proof of Lemma 3.4, so we 
just state the result: 

lim l ° gQL r f' X) =-F(f,\). (3.20) 

As far as the second item on the list is concerned, here we use that by the results of [5] the 
magnitudes of the Gaussian field with distribution (3.18) are positively correlated. (An alternative 
proof of this fact uses reflection positivity.) Invoking the product structure of xa,l and translation 
invariance of Pa, we thus have 

IEa(xa,l)>IPa(|^o| < A) L2 , (3.21) 

where #o is the variable at the origin of the torus. 

It remains to bound Pa(|$o| < A) fr° m below, which we will do by estimating the comple- 
mentary probability from above. We will pass to the Fourier components defined as in the 
proof of Lemma 3.4. Under the measure (3.18), these components have zero mean, the ran- 
dom variables and for different k and k f are uncorrected (a consequence of translation 
invariance), while for the autocorrelation function we get 

^'l^uWi-^ <3 - 22) 

This allows us to use the (quadratic) Chebyshev inequality to derive 

E A (|tfo| 2 ) 1 E A (|^| 2 ) 



> A) < 



A 2 L 2 ^— ' A 2 - /3JA 2 A' 



Inserting this into (3.21) and applying (3.19) and (3.20), the rest of the proof boils down to taking 
logs, dividing by L 2 and letting L — ► oo. □ 

Now we are ready to prove the principal approximation theorem: 

Proof of Theorem 3.1. We just assemble together the previously discussed ingredients. First, our 
constraints (3.7) imply that A < 5 2 and so we can assume that A < ir. Under this condition 
the integrals in (3.3) and (3.13) are over the same set of # r 's an d so by Lemma 3.3 we have the 
uniform bound 

\o g Q LA m__ FLAm < Cl(1+7)/3JA 3. (3.24) 



L 2 

Second, Lemmas 3.4-3.5 ensure that 



nm S up| l0g ^ 2 A((n -F(cf>*)\ < \F(cf>*) - F(<j>*, A)| + log(l - — ^) . (3.25) 
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By the assumptions in (3.7), given an e > we can choose S > such that the right-hand side 
of (3.24) is smaller than e/2. On the other hand, since F(<f>*, A) increases to F(<j>*) as A j and 
since f3A 2 > 1/5, we can certainly choose a A > (satisfying f3JA 2 X > 1) and adjust 5 such 
that also the right-hand side of (3.25) is less than e/2. Combining these observations, the desired 



Remark 1. Physically motivated readers will notice that in both Lemmas 3.4 and 3.5 we have 
introduced a "mass" into the spin-wave spectrum before (or while) removing the indicator xl,a- 
The primary reason for this is the bad behavior of the zero Fourier mode for which the "spin-wave 
Hamiltonian" 1l,4>* provides no decay in the Gaussian weight. 



Having discussed the spin-wave approximations (which will be essential for the arguments in this 
section), we are now ready to start with the proof of phase coexistence. Our basic tool in this 
section will be the chessboard estimates, so we will begin by introducing the notation needed for 
applications of this technique. 

4.1 Chessboard estimates. 

As mentioned previously, in order use chessboard estimates, for technical reasons, we have to 
confine our technical considerations to toroidal geometries. Again we will use to denote the 
torus of L x L sites (as in Section 3 we restrict L to multiples of four). We will consider several 
events which will all take place in a box A# of (B + 1) x (B + 1) sites (which, for definiteness, 
we will assume to be placed with its lower-left corner at the torus "origin"). Since we want to 
be able to cover by translates of A^, we will assume that L is an even multiple of B. Thus, 
if A is an event in A#, then its translate by t\B lattice units in the x-direction and t^B units in 
the y-direction will be denoted by Tt{A), where t = (ti, <2). Here t takes values in a factor torus, 
namely, t £ T L / B . Note that events in the "neighboring" translates of A# may both depend on 
the shared side of the corresponding boxes. 

Let Wl,/3 denote the Gibbs measure on defined from the appropriate torus version of the 
Hamiltonian (1.3) and inverse temperature j3. Specifically, using the "spin-version" of the Hamil- 
tonian (3. 1), the Radon-Nikodym derivative of Wl,/3 with respect to the a priori spin measure Qj L 
is e _/3Wi ( s V-^L,/3' where Z^^ is the corresponding partition function. The statement of the 
chessboard estimates will be considerably easier if we restrict our attention to reflection symmet- 
ric events, which are those A for which S G A implies that the corresponding reflection S* in 
any coordinate plane passing through the center of Ag satisfies S* G A. For these events we will 
also define the constrained partition function 



bound (3.8) is proved. 



□ 



4. Proof of phase coexistence 




(4.1) 
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Theorem 4.1 (Chessboard estimates) Consider the Gibbs measure ^l,/3 ® s defined above. Let 
A\ , . . . , A m be a collection of ( not necessarily distinct) reflection-symmetric events in A b and 
let ti, . . . , t m be distinct vectors from T L / B . Then 



3=1 3=1 



Proof. This is the standard chessboard estimate implied by the reflection positivity condition 
[19-21]. Here we consider reflection positivity in planes "through" sites, which holds in our case 
because we have only nearest and next-nearest neighbor interactions. □ 

Unfortunately, as often happens with chessboard estimates, we may not be able to estimate 
directly the quantity Zi^(A) for the desired event under consideration. Instead, we will de- 
compose A into a collection of more elementary events for which this estimation is easier. Here 
chessboard estimates can be used to establish the following standard (and often implicitly used) 
subadditivity property: 

Lemma 4.2 (Subadditivity) Let the torus Tl and the block Kb be as above and let us consider 
reflection- symmetric events A and (Ak)ke!C m h-B- If A C \J k£l ^Ak, then 

Z L ,M) {B,L? <Y. Z lAM B/L? - (4-3) 

keK. 

Proof. See, e.g., Lemma 6.3 in [3]. □ 

Our succinct recount of the chessboard estimates is now complete. Readers wishing to obtain 
more details on this and related topics are referred to (still succinct) Section 6.1 of [3] or the 
classic references [19-21] and [38]. 

4.2 Good and bad events. 

Here we introduce the notion of good and bad blocks and events. Roughly speaking, a block 
is good if all spins on both sublattices are tolerably close to a Neel state and where the relative 
orientation of the two Neel states is near one of the two optimal values predicted by the spin-wave 
approximation. The bad blocks will of course be all those that are not good. Both these notions 
will involve two parameters: the spin-deviation scale A encountered already in Section 3, and the 
scale k marking the distance to a spin-wave minimum which is still considered good. We will 
keep k small but fixed, while A will have to be decreased (and the block scale B will have to be 
increased, albeit only slowly) as (3 goes to infinity. 

The precise definition is as follows: 

Definition 1 We say that a translate of by Bt, where t £ T^/b, is a good block, or that the 
good block event occurred in this translate if there exist two angles 9* and 4>* such that: 

(1) The angle 4>* satisfies either \4>*\ < n or \(f)* — 180° | < k. 
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(2) The collection of deviation angles # = (i? r ) defined from the angle variables 6 = (9 r ) and 
the angles 0* and <p* via (3.2) obeys 



for all r G T L . 

Let Go be the notation for good-block event with <p* rs 0° and let C/iso be the good-block event 
for </>* in the ^-neighborhood of 180°. The complementary bad-block event will be denoted by B. 
We remark that all these events depend only on the spin configuration (angle variables) in A#. 

Remark 2. It is clear that if either Qq or C/iso occurs (and if k, A <C 1), then the spins in A# are 
indeed well-behaved in the sense of (2.2-2.3) in Theorem 2.1. Explicitly, if r,r' G A# is any pair 
of next-nearest neighbors, then S r ■ S r > is very close to negative one. Moreover, on Go we have 
S r ■ S r > « 1 when r' =r + e x and S r -S r > « — 1 for r' = r + e^, while the opposite relations hold 
on Giso- (Once k, A <C 1, the requisite error is proportional to k for next-nearest neighbors and 
to A 2 for the nearest neighbors.) Thus, the first step in obtaining (2.2-2.3) will be to show that 
any particular block is of a given type of goodness with probability tending to one as f3 —> oo. 

Our goal is to use chessboard estimates to show that, with overwhelming probability, any given 
block is good and that, if one block is good with a known type of goodness, any other given block 
(regardless of the spatial separation) will exhibit the same type of goodness. As it turns out, on the 
basis of Theorem 4.1, both of these will boil down to an efficient estimate of the quantity Z^ p{B) 
defined in Section 4.1. Unfortunately, here we will have to introduce a further partitioning: We 
let Be denote the event that, for some next-nearest neighbor pair r,r' G A#, we have 



This event marks the presence of an energetic "catastrophy" somewhere in the block. As we will 
see, the complementary part of B, 



denotes the situations where the energetics — and the spin-wave approximation — are good but 
where the configuration is not particularly near either of the spin-wave free-energy minima. 

The event Bsw will be further split according to the relative angle between the two near-Neel 
states on even and odd sublattices. Specifically, we let <fr*, i = 1, . . . , s, be s angles uniformly 

(i) 

spaced on the unit circle. Then we let B^ denote the event that the block A^ is bad but such that 
there exists an angle 6* for which the deviation angles # = ($ r ) defined using 6* and 0* = <\>\ 
satisfy \Q r \ < A at each r G A#. (Note that the second part is essentially the definition of the 
good block with the additional stipulation that (p* = in part (1) of Definition 1.) It remains to 

(i) 

show that the B^ indeed cover Bsw- 

Lemma 4.3 Let s be such that sA > Air. Then 



Wr\ < A 



(4.4) 




(4.5) 



B SW = B\ B E 



(4.6) 



s 




(4.7) 



i=i 



ORDER BY DISORDER IN AN 0(2)-SPIN SYSTEM 



17 



Proof. Consider a configuration of angle variables 6 = (9 r ) such that B$w occurs. Since this 
rules out the occurrence of Be, we have 

7r-^<|0 r -0 r ,|<7T + A (4. 8) 

for any next-nearest neighbor pair r,r' G A#. But any two sites on the even sublattice in can 
be reached in less than 73 steps and so 9 T > for any even r' G A# is within A/2 of 9q or #o + tt> 
depending on the parity of r' in the sublattice. Hence, the overall deviations from the appropriate 
Neel state in direction 9* = 9q, where 9q is the variable at the torus "origin," do not exceed A/2 
throughout the even sublattice. Similar considerations apply to the odd sublattice where we use 
the positive x-neighbor of the origin to define the angle 9* + <fi*. 

It remains to show that the above implies that the spin configuration is contained in one of 
the events B<^. Let i = 1, . . . , s be the unique index such that 4>* < < (p* +1 , where <f)* +1 
is to be interpreted as (f>*. Then \<p* — <fi*\ < 2ir/s which by our assumption is less than A/2. 
Consequently, all spins on the even sublattice are within A of either 9* or 9* + tt, depending on 
the parity, while those on the odd sublattice are within A of either 9* + 4>* or 9* + (ft* + tt, again 

(i) 

depending on the parity. In particular, the event B^ occurs, thus proving (4.7). □ 
4.3 Proofs of Theorems 2.1 and 2.2. 

As alluded to in the paragraph before (4.5), the computational part of the proof boils down to 
estimates of the partition functions for events Be and Bsw- These will be provided in next two 
lemmas. We begin with the event Be- 

Lemma 4.4 There exists 5 > and constants C2,cs G (0, oo) such that if 13, J G (0, oo) and 
A G (0,1) satisfy the bounds (3.7), then we have 

lim S up(%^) (B/L)2 < 4B\c 3 (3J) B2 ^e-^ JA2 / B2 . (4.9) 

Proof. When Be occurs, the exists a next-nearest neighbor bond in A# where the associated angle 
variables satisfy (4.5). An easy calculation shows that the energy this bond contributes to the 
Hamiltonian in (3.1) — note that the latter assigns zero energy to the Neel ground states — exceeds 
the J-multiple of 

l + cos(vr- A) =2s in 2 (A). (4 .i ) 

Bounding the sine from below by a linear function, which is justified because A/73 < it, the 
right-hand side is not less than a numerical constant times (A/75) 2 . We thus get 

Z l ^(Be) {B/L)2 < 4B 2 e~ c ^ JA2 / B2 , (4.11) 

where C2 G (0, oo) is a constant and where 473 2 > 273(73 + 1) bounds the number of ways to 
choose the "excited" bond in each translate of A#. 

Our next task is to derive a lower bound on the full partition function. A simple way to get such 
a bound is to insert the indicator that all angle variables 9 T are within A of one of the spin-wave 
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free energy minima, say, 0°. This gives 

^>(^f /2 e- L2F ^° 0) , (4-12) 

where -Fl a is as in (3.3). Fix e > and let S > be as in Theorem 3.1. Then our assumptions 
on (3, A and S and the conclusion (3.5) tell us that 

l T inHZ L ,^ L2 >(^) 1/2 e- F ^. (4.13) 

Let us write the right-hand side as (c3/3J) -1 / 2 , where C3 is a positive constant independent of (3 
and A. Raising this bound to the B 2 power and combining it with (4.1 1) the bound (4.9) is now 
proved. □ 

Next we will attend to a similar estimate for the event Bsw- 

Lemma 4.5 For each k <C 1 and each 7 G (0, 2) there exist numbers p(n) > and 5 > such 
that if A <C k and if (3 J, A and 5 satisfy the bounds in (3.7), then 

limsup ( Z ^} B ^) ) {B/L? < SvrA-V^ 2 . (4.14) 

Proof. Let <p*, i = 1, . . . , s, be s angles uniformly spaced on the unit circle. Suppose that s and A 
satisfy Air < sA < 8tt. In light of the decomposition (4.7) and the subadditivity property from 
Lemma 4.2, it suffices to show that, under the conditions of the lemma, 

^ sup (^f^f^ <e -^ (4 , 5) 

for every i = 1, . . . , s. 

First we note that for 0* nearer than k — A to either 0° or 180° we automatically have C 

Go H Gi80- But then B^ = because the event B^ is a subset of B. By our assumption 
that A <^ k we just need to concentrate only on i = 1, . . . , s such that <fi* is at least, say, k/2 

from 0° or 180°. Here we will use that Z Li/3 (B^ r ) is exactly the {jj) L2/2 multiple of the integral 
in (3.3) with 0* = <j>*, while Zl,p can be bounded from below by a similar quantity for 0* = 
0°,i.e., 

z l ab^)V /l2 



^expj-i^A^ + i^ACO )}. (4.16) 

Let now e > — whose size is to be determined momentarily — and choose 5 > so that 
Theorem 3.1 holds. Then the quantities -Fl a (</>*) on the right-hand side are, asymptotically 
as L — > 00, to within e of the actual spin-wave free energy. Hence, we will have 

limsup ( Zl ^ B ^ V /L < exp{-F(<f>*) +F(0°) +2e}. (4.17) 

L^oo V Z L:j3 ) 

This proves (4.15) with p(n) given as the minimum of F((f>*) — F(0°) — 2e over all relevant i. 
To show that p(n) is positive for k <S 1, we first recall that Theorem 3.2 guarantees that F(4>*) 
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is minimized only by (ft* = 0°, 180°. Since all of the relevant 0* are bounded away from these 
minimizers by at least re/2, choosing e = e(re) > sufficiently small implies p(re) > as 
desired. □ 

Apart from the above estimates, we will need the following simple observation: 
Lemma 4.6 Let re <C 1. Then for any be two neighboring vectors t\,t% G Tl/b> 

r tl (£) nr t2 (g) = (r tl (e )nr t2 (e )) U (r tl (g 180 ) nr t2 (g 180 )). (4.18) 
other words, any two neighboring good blocks are necessarily of the same type of goodness. 

Proof. Since Q = Qq U Gi$o, the set on the right is a subset of the set on the left. The opposite 
inclusion is a simple consequence of the fact that neighboring blocks share a line of sites along 
their boundary. Indeed, suppose the shared part of the boundary is parallel with the y axis. 
For re <S 1, Definition 1 requires that the neighboring boundary spins are nearly aligned in 
a £?i8o-block and nearly antialigned in a C/o-block. Hence, the type of goodness must be the same 
for both blocks. □ 

Now we are ready to prove our main result: 

Proof of Theorem 2.1. As is usual in the arguments based on chessboard estimates, the desired 
Gibbs states will be extracted from the torus measure defined in Section 4.1. Throughout the 
proof we will let (3 be sufficiently large and let A scale as a (negative) power of (3 with exponent 
strictly between 1/3 and 1/2, and B grow slower than any power of (3, e.g., as in 

A = /TT5 and B = \og f3. (4.19) 

We note that these relations (eventually) ensure the validity of the bounds (3.7) for any given 5 > 
and thus make the bounds in Lemmas 4.4-4.5 readily available. 

First we will show that in any typical configuration from most blocks are good. Let t]l 
denote the sum of the ratios on the left-hand side of (4.9) and (4.14), i.e., 

^ Zl,(3 ' ^ Z L,f3 > 

and let r] = limsup^^ t)l- By Theorem 4.1 and Lemma 4.2, the probability of a good block is 
then asymptotically in excess of 1 — rj. On the basis of Lemmas 4.4-4.5, rj is bounded by the sum 
of the right-hand sides of (4.9) and (4.14) which under the assumptions from (4.19) can be made 
as small as desired by increasing (3 appropriately. 

It remains to show that blocks with distinct types of goodness are not likely to occur in one 
configuration. To this end let us first observe that, once re is small, no block can simultaneously 
satisfy both events Qq and (Jiso- Invoking also Lemma 4.6, in any given connected component of 
good blocks the type of goodness is homogeneous throughout the component. (Here the notion 
of connectivity is defined via T L / B , i.e., blocks sharing a line of sites in common, but other def- 
initions would work as well.) We conclude that two blocks exhibiting distinct types of goodness 
must be separated by a closed surface (here *-connected) consisting of bad blocks. 

We will now employ a standard Peierls' estimate. For any t G T L / B the event Qq n T t (Giso) is 
contained in the union of events that the respective blocks are separated by a *-connected surface 
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involving, say, m bad blocks. Using our choice of 77, Lemma 4.2 and Theorem 4.1, the probability 
of any surface of this size is bounded by r\ m . Estimating the number of such surfaces by c m , for 
some sufficiently large c < 00, and noting that m is at least 4, we get 



Obviously, the right-hand side tends to zero as rj J, 0. 

Thus, informally, not only are most blocks good, but most of them are of particular type of 
goodness. To finish the argument, we can condition on a block farthest from the origin to be, 
say, of C/iso-type. This tells us, uniformly in L, that with overwhelming probability the block 
at the origin is of type Giso and similarly for the other type of goodness. The conditional state 
still satisfies the DLR condition for subsets not intersecting the block at the "back" of the torus. 
Taking the limit L — > 00 establishes the existence of two distinct infinite-volume Gibbs states 
which clearly satisfy (2.2-2.3) with e(/3) directly related to rj and the various other parameters (cf 
Remark 2 in Section 4.2). □ 

Proof of Theorem 2.2. This is, of course, just a Mermin- Wagner theorem. Indeed, the Hamil- 
tonian (1.3) satisfies the hypotheses of, e.g., Theorem 1 in [26], which prohibits breaking of any 
(compact) continuous internal symmetry of the model. □ 
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